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Abstract. The decay mode D" — > K^K'^ has no factorizable contribution. We calculate the nonfac- 
torizable chiral loop contributions of order 0{p^) and then we use a heavy-light type chiral quark 
model to calculate nonfactorizable tree level terms, also of order 0{p^), proportional to the gluon 
condensate. Calculated chiral loops and the gluon condensate contributions are of the same order of 
magnitude as the experimental amplitude. 

For nonleptonic decays of D mesons [1] - [10] as well as for ^'s and 5's, the so 
called factorization hypothesis has been commonly used. The factorization hypothesis 
are known to fail badly for nonleptonic K decays [11, 12]. On the other hand, there 
are certain heavy hadron weak decays where factorization might apply. Recently, the 
understanding of factorization for exclusive nonleptonic decays of B mesons in terms 
of QCD in the heavy quark limit has been considerably improved [13]. Following 
[14] we discuss nonfactorizable terms for D decays, in particular for the decay mode 
£)0 ^ ^0^0 In o° K^K^, factorization misses completely, predicting a vanishing 
branching ratio, in contrast with the experimental situation. To see this, note that at tree 
level the — > K^K^ decay might occur due to two annihilation diagrams [1] which 
could potentially create the K'^K'^ state. However, they cancel each other by the GIM 
mechanism. Moreover, in factorization limit, the amplitude is proportional to 

(ifO^>^|0)(0|A''|DO) ^ {pKO-PK^^),fDp'D = . (1) 

In many of the studies (e.g. [2, 3, 4, 5, 7]) this decay has been understood as a result 
of final state interactions (FSI) e.g. [2]. A recent investigation of the K^K^ decay 
mode performed in [3] has focused on the s channel and the t channel one particle 
exchange contributions. 

On the other hand it is well known that factorization does not work in nonleptonic 
K decays. Among many approaches the Chiral Quark Model (xQM) [16] was shown 
to be able to accommodate the intriguing A/ = 1/2 rule in K ^ nn decays, as well as 
CP violating parameters, by systematic involvement of the soft gluon emission forming 
gluon condensates and chiral loops at O (p^) order [12]. In the xQM [16] the light quarks 
(m, J, s) couple to the would-be Goldstone octet mesons (K, 71, rj) in a chiral invariant way. 
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such that all effects are in principle calculable in terms of physical quantities and a few 
model dependent parameters, namely the quark condensate, the gluon condensate and 
the constituent quark mass [12, 16, 17]. 

In the case of D meson decays one has to extend the ideas of the %QM to the sector 
involving a heavy quark (c) using the chiral symmetry of light degrees of freedom as 
well as heavy quark symmetry and Heavy Quark Effective Field Theory (HQEFT). 
Such ideas have already been presented in previous papers [18, 19, 20] and lead to the 
formulation of Heavy-Light Chiral Quark Models (HL^QM). In our formulation of the 
HLxQM Lagrangian, an unknown coupling constant appears in the term that couples 
the heavy meson to a heavy and a light quark [14]. Our strategy is to relate expressions 
involving this coupling to physical quantities, as it is done within the XQM [12]. We 
perform the bosonization by integrating out the light and heavy quarks and obtain a 
heavy quark symmetric chiral Lagrangian involving light and heavy mesons [21, 22, 23]. 



Because the Oip) (factorizable) contribution is zero as seen in Eq. (1), we approach 
to the K^K^ decay by calculating systematically O (p^) contributions. We do this 
by including first the nonfactorizable contributions coming from the chiral loops. These 
are based on the weak Lagrangian corresponding to the factorizable o{p) terms for 
Z)0 ^ n+K- andZ)0 ^ K+K . Second, we consider the gluon condensate contributions, 
also of O (p^) within the xQM and HL%QM framework. The energy release in D — > KK 
is p = 788 MeV and hence p/A^ (for A,^ > 1 GeV), is close to unity. The next to leading 
0{p^) terms might be almost of the same order of magnitude compared to our 0{p^) 
terms. However, we expect a weak suppression of the order p^/A^. On the other hand, 

the inclusion of O (p^) order in this framework is not straightforward. Before doing loop 
calculations at that order, one has to find a reliable framework to include light resonances 
like p, K*, ao(980), /o(975) etc. The poorly known scalar resonances would introduce a 
rather large uncertainty [3]. Right now, a consistent calculation of this or higher orders 
does not seem to be possible. 

Note that we have omitted 1/mg terms in the framework of HQEFT. 

The effective weak Lagrangian at quark level relevant for D — > %%,KK is 



are quark operators. 

Using Fierz transformations [14] one obtains operators Qa = Qg^/Nc+R 



Qc = Qb /Nc + Rf, Qf = Qa/Nc + Ra and 27 = Qc/Nc + Rc, where the 



R's correspond to color exchange between two currents_ and are genuinely 

nonfactorizable: Ra = 2{sLY't" ci^){uLt''y^SL), Rc = 2(dLtt''cL){uLt''y^dL), 
Rb{i) = "^iui^t" CL){^i^t"y^qL) , {q = s,d). The operators can be written in 
terms of currents [14]. The factorization approach amounts to writing the currents in 




Lw = G ca{Qa-Qc) +cb{Q, 




(2) 



where G = -2\/2GfK^V;^, and 

Qa = (sLtcLl {uanSl) , Qc = {dLfci) {uL^dl) , 

= (uLtcL) {qanqL) , {q = s,d), 



(3) 
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FIGURE 1. The diagrams which give nonzero amplitudes. 



terms of hadron (in our case meson) fields only, so that the operator 2^ — ^he 
left equation is equal to the product of two meson currents. The color currents are then 
zero if hadronized (mesons are color singlet objects). There is also a replacement of the 
Wilson coefficients in the hadronized effective weak Lagrangian ca.b ca,s(1 + 1 /Nc). 
Combining heavy quark symmetry and chiral symmetry of the light sector, we can 
obtain the weak chiral Lagrangian for nonleptonic D meson decays due to factorizable 
terms. Then we can first use this to calculate nonfactorizable contributions due to 
chiral loops. Second, we can calculate the color currents' contribution using the gluon 
condensate within the framework of the HL%QM. 

Treating the light pseudoscalar mesons as pseudo-Goldstone bosons one obtains the 
usual 0{p^) chiral Lagrangian [14] and from this Lagrangian, we can deduce the light 
weak current to o{p). 

In the heavy meson sector interacting with light mesons we have used the lowest order 
0{p) chiral Lagrangian (for details see [14]). Using symmetry arguments, the heavy- 
light weak current is bosonized to O (p^), with the unknown coupling an related to the 
physical decay constant fo- 

In the factorization limit there are no contributions to K^K^ at tree level. The 

observation of a partial decay width B{D K^K^) = (6.5 ± 1.8) x 10""^ on the other 
hand implies that we can expect sizable contributions at the one loop level. Calculations 
to one loop in the framework of combined chiral perturbation theory and HQEFT 
involves a construction of the most general effective Lagrangian that has the correct 
symmetry properties in order to make the renormalization work. Our calculations were 
done in the strict MS renormalization scheme. 

Writing down the most general one loop graphs with two outgoing Goldstone bosons 
(K'^ and K^) one arrives at 26 Feynman diagrams. The expressions for nonzero ampli- 
tudes corresponding to the graphs on Fig. 1 are given in [14]. 



FIGURE 2. Feynman diagram for bosonization of left-handed current to order o {p) 
The partial decay width for the decay — > K^K^ is then 

where F = J^n^n is the sum of the ampUtudes corresponding to the graphs on Fig.l and 
p is the K'^{K^) three - momentum in the rest frame. 

In numerical calculation we use the values of an, g and / (/ is related to the % 
meson decay constant) obtained within the same framework [21, 23, 27, 28, 29]. The 
coupling g is extracted from existing experimental data on D* Dn and D* Dy 
decays. The analysis in [23] includes chiral corrections at one loop order and yields 
g — 0-27^Q Q2-o 02' leaving the sign undetermined. Recently CLEO Collaboration has 
measured the Z)*+ decay width [25]. By combining this result with existing data on D* 
decay widths [15], we obtain value g = 0.57 ± 0.08. We present results for g — 0.27 and 
g = 0.57. The larger value seems to be in better agreement with the results coming from 
different approaches listed in [21]. We put everywhere = 1 GeV ~ Aj^. 

For the Wilson coefficients ca,b of (2) we use ca = 1 . 10 ± 0.05 and cb = —0.06 ±0.12 
[24], calculated at the scale // = 1 GeV with the number of colors Nc = 3. Due to the 
suppression of cb in comparison with ca, we do not include terms proportional to cb- 

TABLE 1. Table of the one chiral loop amphtudes (see Fig. 1), 
where M =Y.n^n- The second column shows the ampUtudes cal- 
culated using g = 0.27 while the third column amphtudes have been 
calculated using g = 0.57. 



!Mi[x 10-^GeV] {g = 0.27) Mi[x 10-^GeV] {g = 0.57) 



Ml 


-0.42 


-0.82 


Ml 


-0.31 


-0.62 


M'i 


-0.62 


-1.23 


M4 


0.75 -2.54 i 


0.70 -2.37 i 


Ms 


-0.81 


-0.76 


Me 


-0.61 


-0.57 


M-] 


-0.99 


-0.92 


Ms 


0.91 


0.85 




-2.11 -2.54 i 


-3.37 -2.37 i 



The imaginary part of the amplitude comes from the F4 graph, when the TC's or the 
^'s in the loops are on-shell. All other graphs contribute only to the real part of the 
amplitude. The imaginary part of the amplitude is scale and scheme independent within 
chiral perturbation theory. This amplitude is also obtained from unitarity, and is valid 
beyond the chiral loop expansion. 
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FIGURE 3. Diagram for bosonization of the color current to O (p^) 
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FIGURE 4. Diagrams representing bosonization of heavy-light weak current. The boldface line repre- 
sents the heavy quark, the sohd hne the hght quark. 



In the effective weak Lagrangian there are, after Fierz transformations, terms that 
involve color currents. As mesons are color singlet objects, the product of color currents 
does not contribute at meson level in the factorization limit. However, at quark level 
they do contribute through the gluon condensate. In order to estimate this contribution 
we have to establish the connection between the underlying quark-gluon dynamics and 
the meson level picture. This is done through the use of a Heavy-Light Chiral Quark 
Model (HL%QM). Our starting point is the Lagrangian containing both quark and meson 
fields [14, 18, 19, 20, 26]. After bosonizing the light weak current at the order O (p) and 
0{p^) (see Fig. 2 and Fig. 3) as in [14] one should bosonize the heavy - light weak 
current integrating out quark fields as presented on Fig. 4. The product of two external 
gluon fields (G in Fig. 3 and Fig. 4 ) is interpreted as the gluon condensate ((G^)) [14]. 
This contribution is of the 0{p^) order. In the language of chiral perturbation theory, 
the divergent part of the counterterm has the Lorentz and flavor structure of the effective 
Lagrangian given in eq. (47) in [14]. 

By taking into account various relations of the loop integrals we determine 



which is of the same order of magnitude as the chiral loop contributions in Table 1. 
Adding both the chiral loops and the gluon condensate (5) contributions, we obtain the 
total amplitude to O (p^) 



M (Z)° ^ ^^i^0)(G2) ^ 0.43 X lO-'^GeV ; 



(5) 



g = 0.27; iM'Th= (-1.68 -2.54/) x lO-'^GeV 
g = 0.57; fAfxh = (-2.94 - 2.37/) x lO-'^GeV . 



(6) 
(7) 



or in terms of branching ratio 



^ = 0.27; 5(DO^i^o^O)Th= (4.2±1.4) x 10"^ 
g = 0.57; B{D^ f^K%t, = (6.5 ± 1.7) x lO'^ 



(8) 
(9) 



where the estimated uncertainties reflect the uncertainties in the rest of input parameters. 
These results should be compared with experimental data [15] 5(Z)^ — > K'^K^) — (6.5 ± 
1.8) xlO-'^. 

We can summarize that we indicate the leading nonfactorizable contributions to — > 
K^K'^. Even though the use of chiral perturbation theory in this decay mode could 
be questioned, the calculated chiral loops can be considered as part of the final state 
interactions. In the treatment of the final state interactions the light pseudoscalar meson 
exchanges have to be present due to unitarity. Although, the next to leading o{p^) 
order terms might give sizable contributions to this decay, we have demonstrated that 
contributions due to the chiral loops and gluon condensates are of the same order of 
magnitude as the amplitude extracted from the experimental result. 
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